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ABSTRACT 


The  effect  of  the  wall  on  the  motion  of  elongated  aerosol 
particles  was  calculated  on  the  basis  of  the(powerful)  slender 
body  theory  for  a  quiet  and  flowing  medium.  The  cases  treated 
were  those  of  cylindrical  particles  moving  in  a  quiet  environ¬ 
ment  or  a  Couette  or  a  Poiseuille  flow  near  plannar  surfaces  with 
various  inclinations  to  the  horizontal. 

Experiments  of  trajectory  photographing  were  performed  in  a  liquid- 
tank  under  creeping  flow  conditions  and  typical  (aerosol)  Reynolds 
numbers  whereby  very  good  correspondence  with  computations  was 
obtained . 

The  Brownian  diffusion  of  (very  small)  cylindrical  or  discoidal 
aerosol  particles  was  also  studied  while  being  based  on  Brenners' 
formal  structure.  To  evaluate  the  signficant  orientation  density 
function,  the  (mid-diameter) rotational  diffusion  coefficient  of  the 
particles  was  calculated  by  a  statistical-mechanical  method.  The 
values  of  the  coefficients  were  found  to  be  larger  than  those  computed 
on  the  basis  of  continuum  theory  by  a  factor  of  about  30. 

The  typical  diffusion  situation  of  deposition  on  a  plane  within  a  quiet 
ai r  was  treated. 
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INTRODUCTION 


During  the  past  several  years  a  series  of  studies  on  the  transport 
mechanics  of  non-spher ical  aerosol  particles  was  carried  out  mostly 
with  the  support  of  the  European  Research  Office,  U.S.  Army. These  studies 
were  summarized  in  the  final  reports: 

The  Dynamics  of  Nonspherical  Particles 

I.  Controlled  Generation  of  Nonspherical  Aerosols  and  Methods 
for  Their  Size,  Concentration  and  Electric  Size  Measurement  (Contract 
DAJA  37-72-3875;  1974). 

II.  Mobility  in  Non-Orienting  Fields  (Contract  DAJA  37-74-C- 

1208;  1975). 

III.  Mobility  and  Deposition  in  Still  Air  (Contract  DAER0-75-G-021 ; 
1977) • 

IV.  Deposition  from  Flow  (Contract  DACR0-76-G-014;  1978) 
and  in  the  parallel  journal  publications: 

1.  "On  the  stochastic  Motion  of  Nonspherical  Aerosol  Particles. I .  The 
Aerodynamic  Radius  Concept"  by  Isaiah  Gallily  and  Aaron  Hi-Cohen,  J.  Colloid 
Interface  Scl.  5j6,  443-59  (1976) 

2.  " _ , '  I.  The  Overall  Drift  Angle  in  Sedimentation",  by 

Matitia  A.  Weiss,  Aaron-Hi  Cohen  and  Isaiah  Gallily,  J.  Aerosol  Sci.  9> 
527-41  (1978). 

3.  "On  the  Orderly  Nature  of  the  Motion  of  Nonspherical  Aerosol 
Particles.  I.  Deposition  from  a  Laminar  Flow."  by  Isaiah  Gallily  and 
Alferd  D.  Eisner,  J.  Colloid  interface  Sci.68,  326-37  (1979). 
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Likewise,  they  were  reported  in  the  following  articles: 

A,  "On  the  Orderly  Nature  of  the  Motion  of  Nonspherical  Aerosol 
Particles.  II.  Inertial  Collision  between  a  Spherical  Large  Drop¬ 
let  and  an  Axially  Symmetrical  Elongated  Particle"  by  Isaiah  Gallily 
and  Aaron-Hi  Cohe,  J.  Colloid  Interface  Sci.  68^,  338-56  ( 1 979) - 

5.  " _ ,11. _ ,  Further  Results,"  by  Isaiah  Gallily  and 

Aaron  Hi-Cohen,  ibid.  7J_,  628*30  (1979). 

The  efforts  invested  in  all  of  these  investigations  culminated  in  the 
establishment  of  a  method  of  claculation  of  the  gravitational  sedimentation 
in  still  and  laminarily  flowing  gas,  and  of  the  inertial  impaction  on 

various  surfaces  of  ax i -symmetrical  particles  such  as  cylinders  and  i 

disks.  The  method  has  been  ellaborated  in  the  (above)  publications  1-5  and 
partially  checked  by  experiments;  however,  it  is  far  from  being  complete. 

One  problem  in  the  analysis  of  our  particles'  transfer  phenomena  which 
has  not  been  considered  is  the  effect  of  the  surface  exposed  to  the  aerosol  on 
the  collection  process  itself.  The  effect,  related  here  only  to  the  motion 
of  the  particles  in  the  adjoining  gas  layer,  is  due  to  . 

a.  The  f 1 u i ddynami c  interaction  between  a  particle  approaching  a 

surface  and  the  latter  which  essentially  stems  from  the  change  of  the  flow  field 
due  to  the  emergence  of  an  additional  boundary  condition. 

b.  The  presence  of  long  (electrostatic)  or  short  (Van  der  Waals)  surface 
particle  forces. 

Another  untreated  problem  concerns  the  stochastic  diffusional  turbulent  and/or 
Brownian  movement  of  the  particles. 


A 
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It  is  needless  to  stress  the  significance  of  the  first  problem  which  has 
been  treated  scarcely  for  the  case  of  spherical  particles  by  Brenner(l), 
Gal  lily  and  Mahrer  (2)  and  others  and  has  been  almost  untouched  in  the 
case  of  the  nonepherical  ones. 

In  most  studies  on  the  orderly  motion  of  any  particle,  expressed  by  the 
equation  of  translation 


™  Au./ dr  =  Z  Fe,.  +  F, 

t 


M 


and  rotation 


d(l"-*>)/dt  =Z  Lt+  Lf  [2] 

the  fluiddynamic  operating  force  Fp  and  torque  Lp  have  been  calculated 
'or  the  substitute  situation  in  which  the  particle  is  infinitely  far  from 
the  (col lect ing) surface.  This  has  been  obviously  a  simplification  whose 
justification  was  checked  in  the  calculation  of  Fp  and  Lp  for  some  spheri¬ 
cal  cases  (1,2).  For  the  nonspherical  case,  it  has  not  been  studied  at  all. 

The  importance  of  the  second  problem  is  self  evident  for  sufficiently 
turbulent  field  or  very  small  particles  for  which  the  fluctuating  part  of  the 
force  Fp,  viz.  Fp  -  (fp)  *  can  not  be  neglected  in  comparison  with 


the  mean  one  <F_>. 


AIM  OF  STUDY 


In  view  of  the  above  review  of  the  state  of  research,  it  was  thought 
worthwhile  to  investigate: 

1.  The  effect  of  the  collecting  surface  on  the  motion  of  nonsphe- 
rical  particles  near,  and  their  deposition  on  it,  and,  independently, 

2.  The  diffusion  of  nonspherical  particles. 

These  two  problems  seemed  immediately  too  formidable  for  a  general 
treatment;  thus  their  study  was  restricted  in  the  first  subject  to  the 
case  of  cylindrical  particles  much  larger  than  the  mean  free  path  1  of 
the  medium's  molecules  (  K-m  <^<  1  ),  and  a  planar  solid  wall,  and 

in  the  second  subject  to  the  Brownian  diffusion  of  cylindrical  and  discoidal 
particles  much  smaller  than  L  (  K<*  <<  1  )  .  Cylindrical  and 

discoidal  forms  respectively  serve  as  idealized  model  shapes  of  the 
real  (harmful ) asbestos  or  other  significant  fibrous  particles  and  of 
latelet-like  ones.  So,  the  restriction  to  the  above  morphology  is  not 
very  limiting  in  significance. 


A.  THE  WALL  EFFECT 


I.  Theoretical 
1 .  Method 

a.  General 


Here  we  preferred  to  use  slender  body  theory  rather  than  the  (limited) 
method  of  Wakiya  which  is  based  on  Fourier  analysis  technique  (3). 


The  slender  body  theory  method  is  much  more  general  than  the  latter; 
it  can  be  employed  in  situations  where  the  (aerosol)  particle  is  even 
farfrom  walls  with  a  better  advantage  than  that  obtained  by  being  based 
on  the  much  used  Oberbeck's  (A)  and  Jeffery's  (5)  equations.  In  other 
words,  it  is  very  useful. 

The  method  is  apt  for  the  treatment  of  the  motion  of  particles  whose 
characterist  ic  size  rc  (=  V  £  /  ^  tt  )  is  much  larger 

i 

than  the  mean  free  path  1  of  the  mediums  molecules 

and  whose  characteristic  Reynolds  number  Re  is  much  smaller  than 

unity.  In  fhat  case,  the  equations  of  fluid  motion  to  be  solved  are 


p-  y1  v  ("»■)  -  7  f>  r  0 


[31 


and 

7«V(r)  s  0  .  ^ 

The  theory,  developed  by  Batchelor  (6)  ,Chwang  and  Wu(7)  ,  Blake  (8)  , 

Cox  (9)  and  others,  deals  with  elongated  particles  in  which  the  ratio  k 
between  the  characteristic  cross-sectional  radius  b  and  the  semi-length  1 
is  W  =  1  ,the  change  of  the  local  equivalent  cross-sectional 

*  For  symbols  see  NOMENCLATURE 


rad i us 


is  small  compared  with  &>>  /  (r  »  and  the  local  radius 


Xi 


Fig.  1.  A  slender  body  (not  necessarily  a  straight  cylinder); 
coordinates  used  in  the  analysis. 


Under  these  conditions,  the  disturbance  velocity  v  (r)  produced  by  the 
moving  particle  is  given  by  the  multi  pole  expansion 


7 


t  f 

(O  =  \  ij  u  I(,  „  r,  r']  dr1  +  j"  /jit  (r1)  Uijb  (r, 4  r‘ 

-t  -( 

+  -fj.*  (*-')  U  |»  *  (  r,  y~')  clr‘  +  .  .  -  ^ 


where,  with  the  Einstein's  convention  adhered  to,  ^a(r')  (Stokeslet)i 


••  ^ r'  ^ 


(doublet)  , 


(  r1  ) 


etc.  are 


(imaginary)  poles  causing  the  disturbances  u  j  (r,  r‘)  etc,;  the  latter 
are  actually  the  proper  Green  functions  (APPENDIX  1). 


The  calculation  of  the  (critically  required  )  fluiddynamic  force  F p 
and  torque  L  p  appearing  in  Eqs.  [1]  and  [2]  can  be  performed 
once  the  Stokelet  -f  (r')  distribution  is  known  since  (Fig.  1 ) 


ff  =  (  fcs)  d* 


cF  =  J  ( X  x  4  (s))  d  s 

-l 


The  finding  of 


(r‘)  cannot  be  carried  out  by  the  solution  of  the 


full  Eq.  [5]  with  the  proper  boundary  conditions  on  the  particle 
(v  (r)»  U  (r)  ■t  U  (r)  +  cu  x  x  for  r  — »  S  +  b)  and 

on  the  wall  (v(»)  =  0  )  •  Thus,  one  uses  as  an  approximation  for  t  »  4 

(Fig.l)  an  "outer  expansion"  in  which 


8 


l 

V;  fr)  =  j\:  '  rr‘  u.  ij  (  r,  r'  0/  *•'  C8l 

-l 

only,  and  for  j-  <«.  J  1  (for  which  the  disturbance  can  be  assumed  to 
evolve  from  a  portion  of  an  infinite  straight  cylinder  of  radius  <J-r, 
an  "inner  expansion"  (6) 


v, 

1 


1 


2r r  u 


c;  O)  [f*j 


i, 

§ 


!<(s)]  4  Uy(<) 


[91 


and 


vt.  Cr) 


_ 1 _ 

^TTp. 


C' (r)[6*i 


-^-4  +  [10] 

$ 


where  I*  ,  I’  (  =  2  ,  5  ) 


respectively  stand  for  the  particle- locked 


local  coordinates  x,.  ,  xx>  ,  x3i  ,  C(  ,  C  '■  are  quantities 

J 

affected  by  the  wall  and  K  ,  K^;  are  constants  defined  by  the 

a 

geometry  of  the  cross-section  at  e'  (Fig.  1).  This  "inner  expansion" 
can  be  also  written  in  the  vectorial  form  fit  for  any  coordinate  system 


=  M-C  &i-|[*+M.Tr.C4l-C-I*C  +  ^n-f<-U  [tl] 

where  H  ,  C,  "IT,  I,  I'  are  defined  in  APPENDIX  1. 

The  two  approximations  are  matched  In  a  region  in  which  they  are  both 
assumed  to  hold  ( &  «  $«  £  )  ,  and  the  relevant  constants  so  derived. 


In  the  case  of  the  "outer  expansion",  one  can  write 


...  Si--'.:: 


(  r ^  r*  )  =  U  ,j  C  ^  »  r  / 


where  u.,*”,  (  r  ,  r' j  has  singularities  on  the  particles  axis 

and  -*■  C  r?  r*  .  is  non-singular  when  §  » 


b.  Straight  Bods 


In  this  case  Eq.  [8]  i  s  integrated  by  Batchelor's  approximation  to  ['ll]  (&) 


which  gives 


kitfi-y/  (r)  =  {~r  +  [- 


-I  V  +  —  M-  f?  +^V.  4  J-«f  [13] 

l",R  and  _A.  being  defined  in  APPENDIX  1  (  <T  =  [  (2  t/(r)l  ' 

and  X  a  (&»>•/  6 r  )  .  Now,  matching  the  coefficients  of  the 

(singular)  part  in  [11]  and  [13]  one  gets 

§  ’ 

C(.)  =  Us)  ,  C14J 


and,  doing  the  same  for  the  non-singular  part, 
hi t|lx  U($)  •»  ^  1  "  ^  "f  ■*  [15] 

where  l‘=-l"+4  l'  and  TT  0  -  7T  •  M 
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At  this  stage  Ci)  can  be  numerically  calculated  from  [15] 
by  dividing  our  particles'  axis  into  N  equal  parts,  to  each  of  which  is 
associated  an  average  ;  (s)  ,  and  performing  the  integration  of 

aVi  (  r  i  r  }  and  respectively  appearing  in 

JV,  and  R  ;  so,  for  the  two-dimensional  case  of  particles'  motion, 

one  obtains  2  N  equations  for  2  N  unknowns  while  for  the  three- 
dimensional  case  one  obtains  3  N  equations  for  3  N  unknowns. 

According  to  Batchelor  (6),  it  is  possible  to  use  here  a  quick  iterative 
method  by  writing 

4.  CO  =  -►  +  (fa-frpU)  +  ••• 

i  t  *  '»  • 

1  -Ibe.WaliJffect 

Dealing  with  a  cylindrical  particle  undergoing  a  translational 
and  rotational  motion  near  a  planar  solid  wall  (Fig.  2)  and  statisfying 


Fig. 2.  A  slender  cylinder-planar  solid  wall  system,  coordinates 
used  in  analysis  , 


the  no-slip  condition  on  all  boundaries,  one  has  to  introduce  on  the 
reflected  Image  of  the  particles'  axis  a  line-distribution  of  point 
forces,  -.f,(x)  (a  line  distribution  of  doublets,  » 

and  a  line  distribution  of  quadruples,  -  (x)  (8) 
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In  this  case 


<»  r  Sit 


If  -  r<  |* 


D71 


and 


_  1  N  _ 

Scj 

(r;  - 

r.')(rv- rf  ') 

.  C  v 

’ 

r  ;  =: 

|r*i 

3k 

n*  r; 

&i.  r?  _ 

!r*l* 

1 

r*lr 

_ 

ir*l 5 

oc  i  s 

1 

or  2 

and 

L  i  « 

/  4  • 

3  n-  •  r„  •  r, 


.•I*- 


1,2,3  respectively 


ii.  Mot ion_wi thin_a_Flow 

Similarly  to  the  "outer  expansion"  In  a  quiet  medium.  It  Is  possible 
to  write  in  this  case  for  &/$  ^  ^ 


V 

v'Cr)  =  |  ^  O)  u y  (r,r‘)  of**' 


P9] 


■l 


in  which  V(r)sVCr)-  V*(r)  and  V^Cr)  Is  the 

undisturbed  flow  ((v(r)  — ►  v\r)  cs  r - ♦  )  ;  the 

"Inner  expansion"  can  be  expressed  by  an  equation  similar  to  [ifl 

but  where  C  =  (  C  ,  C'a  ,  CJ  )  ,s  now  affected  b*  the 
flow. 


-  >3  - 

Thus,  after  performing  the  necessary  matching  between  [19] and  the 
modified  [11’,  one  gets  like  [l5[ 

■rr  k  U'  (s)  a  C -jp  +  ?-r  [-  ~  ~ ■■ '  ■■  ]^H  ■  4  (»)  +  I<9>.  {  O 

-  TT  •  M  •  -f  -  si  ■+  -  M  •  R  +  -A.  (S[  [20] 

in  which  U'Cs'/  =  U  (*  /  -  ,  U  Is  defined  by 

v(r) - }  U  Cf)  ^or  r  - »  S  b  ,and  all  other  symbols  are 


as  marked  above. 


Expanding  v*0)  on  the  (straight)  particles'  axis  in  the 


+  ,  N  -CS  («)  Q  O  _  OO  * 

v.  (x)  =  4-  X  ■+  JB i  X 


and  putting  xK=xi>t  (i^  being  oos  6  and  iz  being  sind  ), 
one  has  for  U'  Cs) 

Ui  (X)  =  U,(x)-v7(x)  =  V.  -  B;X  ...  [22] 


where  £.  =  €  :>  u  cj«  i w  -  .£  ■  .  W  :  Is  the  rotational  velocity 

t>  Q  0  W  1  J 

of  the  particle  and  V-  is  its  translational  velocity. 
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2.  Computational  Procedure 

The  trajections  of  our  particle  were  derived  from  the  solution  of  Eqs.  .1. 
and  '2'.  on  the  specific  assumptions  that  the  particles  are: 

i.  homogeneously  dense 

and 

ii.  affected  only  by  the  gravitational  force  mg,  the  fluiddynamic 

force  Fp  (Eq.  *6]  )  and  the  fluiddynamic  torque  Lp. 

At  first  were  treated  two-dimensional  cases,  In  which  particles  axis 
stays  in  a  plane  containing  g  ;  then  we  considered  three  dimensional 
si  tuations. 

For  the  two-dimensional  cases,  where  the  [1]  and  [2]  turn  into 


yr>  d  V//  dt  =  (  -  -m' )  3  +  F  f 

and  (since  go,  =  coz  -  o; 

co5  =  L5  , 


[23J 


[24J 


we  employed: 

i.  A  quick  technique  based  on  Eq.  [16] 

and 

ii  An  elaborate  technique  described  in  APPENDIX  1. 

The  computation  itself  was  performed  {on  the  Hebrew  University  CYBER  7**  CDC  com 
puter)  by  the  use  of  a  fourth-order  Range  -Kutta  method;  the  time  intervals 
were  from  5  •  10_it  to  3  •  10_i*  second,  depending  on  the  (decreasing) 


gap  between  the  particle  and  the  wall. 
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II.  Experi cental 

The  experiments  were  essentially  intdnded  to  check  the  calculations  in 
their  significant  points.  Thus,  one  could  in  principle  either: 

1.  Determine  the  force  Fy  and  torque  Lp  from  measured  trajectories  and 
orientations  and  compare  with  theory 

or 

2.  Determine  the  whole  trajectories  of  the  particles  with  their  orienta¬ 
tions  and  test  against  computations. 

At  any  case  it  was  necessary  to  measure  particles'  location  and 
orientation  as  a  function  of  time;  so,  since  the  observation  of  small  cylind¬ 
rical  aerosol  particles  is  complicated  by  light  scattering  phenomena 
(10),  it  was  thought  worthwhile  to  use  large  (macro)  particles  which  move  in  o 
in  a  creeping  motion  {  Re  <<  l  }  too. 

1 .  Method 

In  essence  and  without  loss  of  generality,  our  experiments  were  carried 
photographing 

out  by  /a  two-dimensional  motion  of  a  homogeneously  dense  cylindrical  body 
falling  under  the  action  of  gravitation  in  an  oil  medium.  Here,  out  of  the 
above  two  alternatives,  we  chose  the  latter  which  Is  more  exact  than  the 
former  as  it  does  not  involve  differentiation  of  location  and  velocity  data 
curves. 

2.  Apparatus 

The  apparatus  designed  and  constructed  for  the  above  purpose  (Fig.  3, a 
and  b)  consisted  of: 


a*- 
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A  tank  A  containing  silicone  oil  in  which  was  suspended  a  polished  steel 
plate  B,  a  device  C  to  release  a  metallic  cylinder  C  into  the  liquid,  a 
movie-camera  E  and  light  source  for  the  trajectory  photographing,  a  moving 
platform  F  for  the  camera,  and  a  scale  G  -mirror  -  light  source  combination 
for  locating  the  particle. 

The  tank  was  made  of  Plexiglas  in  the  (inner)  dimensions  of  50  cm  height, 

50  cm.  breadth  and  3  cm.  width  and  was  kept  vertical  by  Its  own  weight. 

The  liquid  (F  111/10000  dimethyl  silicone  fluid,  ICI,  England)  had  a 

L 

(nomi nal ) vi scos i ty  of  10  centistokes  or  about  100  poise  dynamic 

viscosity  at  25°C  with  a  temperature  coefficient  of  0.02V,,/oC,  and  a 
density  of  O.968  g.  cm  ^  at  25°C.  This  li  quid  was  satisfactorily  Newtonian 
as  tested  before  each  experiments  by  measuring  p  at  various  shear  rates 
3  v i  /  3  xj  by  the  settling  velocity  of  various  sized  Stokesian  steel 
balls.  The  cylinder  used  as  a  tested  particle  had  a  radius  b  of  0.0075  cm, 
a  semi-length  1  of  1.5  cm.  and  a  density  of  7-780  g.  cm  ^ .  It  was  intro¬ 
duced  into  the  liquid  by  a  slow  sliding  movement  in  the  groove  of  the  devices' 
edge  K  (Fig.  3  b)  and  along  the  middle  (vertical)  line  of  the  plate 
(wall).  The  camera  was  a  Bolex  H  16  Reflex  movie  instrument  (Baillard, 
Switzerland)  equipped  with  a  1:  18  lens. 


V 


■.  y  is. . 
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3.  Procedure 

The  viscosity  jj-  of  the  liquid  was  determined  before  every  run  as 
described  above;  it  was  found  to  be  constant,  within  the  experimental 
error,  at  the  value  of  1 1 2  +  2  poise  for  16.4°C.  The  temperature  of 
the  oil,  measured  at  various  locations  in  the  tank,  was  constant  to  a 
degree  of  0.02°C,  which  amounted  to  a  viscosity  change  of  0.02  poise. 

The  trajectory  photographing  was  performed  while  keeping  the  center  of 
the  falling  cylinder  on  the  optical  axis  of  the  camera  as  well  as 
possible.  The  reading  of  position  and  orientation  of  the  particles  image  on 
the  movie  film  was  carried  out  with  a  Nikon  6CT2  profile  projector  (Nikon, 
Japan)  to  a  (nominal)  accuracy  of  +  0.001  mm. 

Here  it  should  be  pointed  out  about  the  critical  requirement  of  two- 
dimensionality  to  which  a  self-restriction  was  imposed.  Due  to  the  cylindri¬ 
cal  rod  release-operation, there  was  always  a  chance  that  the  particle  will 
not  move  with  its  axis  in  the  plane  of  g  .  So,  in  every  experiment,  the 
film  image  of  the  falling  rod  was  aposteriori  checked  against  the  expected 
value  deduced  from  the  actual  size  and  orientation  of  the  body,  and, when 
a  deviation  occured,this  run  was  discarded. 

III.  Results 

As  the  checking  experiments  had  to  be  restricted  to  a  liquid  medium 
and  macro  particles,  the  calculations  were  confined  to  that  situation  too. 
Dealing  with  physical  cases  which  can  be  analyzed  by  continuous,  classical 
fluiddynamic  methods,  this  did  not  matter  as  long  as  creeping  motion 
conditions  (  Re  < <  1  K-n  <<  1  )  still  prevailed  and  the  liquid 


specific  effects  associated  with  the  acceleration  of  the  particles,  viz.  the 


19 


added  mass  and  moments  of  inertia  and  the  Basset  term  influence  (11) 
could  be  neglected. 

The  cases  selected  for  computation  were  chosen  to  be  a  two-dimensional 
motion  in  a  still  environment  near: 

i.  A  vertical  wall  ^cigs.  5  through  8), 

ii.  An  horizontal  wall  (Figs.  9  through  11), 

iii.  An  inclined  wall  (Figs.  12  through  19) 

and  a  two-dimensional  motion  within: 

i.  A  two-dimensional  Couette  flow  over  an  horizontal  wall  (Fig.  20) 

ii.  A  two-dimensional  Poi6euilie  flow  through  an  horizontal  channel 
(Fig.  21). 

All  of  these  cases  simulate  real  situations.  Thus,  for  example,  if  the 
particles  are  small  enough,  then  any  corrugated  surface  can  be  represented 
by  inter-connected  plan  ar  portions  of  various  inclinations  to  the  horizontal 
likewise,  a  general  laminar  field  can  be  approximated  by  a  series  of  inter- 
adjoining  Couette  flows.  The  restriction  to  two-dimensional  motion  which 
is  relaxed  latter  on  in  a  further  study  (not  reported  here)  is  not  serious 
since  in  many  instances  such  a  situation  is  indicative  enough  of  the  physical 
occurences . 

The  distances  and  angles  employed  in  the  calculations  are  brought  out 
in  Fig.  k  (the  unit  of  distance  is  always  the  semi-length  1).  Thus  the 
undisturbed  Couette  flow  was  expressed  by 


(x,)  =  k,  I?, 


[25] 
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where  «•  is  the  (constant)  velocity  gradient,  while  the  (undisturbed) 
Poiseu.ille  flow  was  given  by 


( x ,  . 


26: 


where  is  a  (flow)  constant  and  a  is  the  half  width  of  the  channel. 

The  typical  results  shown  here  are  presented  in  terms  of  the  (operationa¬ 
lly)  important  drift  d i splacements  and  linear  velocities  and  the 
intuition-assisting  rotational  velocities  and  particle  orientation. 


Fig.  4.  A  slender  body  near  a  wall.  Distances 
and  angles  used  in  calculations. 


xig  i /sec) 


ig.  9.  Vertical  velocity  of  a  cylinder  falling  towards 
an  horizontal  plane  as  a  function  of  Dj  and  for 
various  initial  orientations  0O  (in  rads). 

Same  cylinder's  and  liquid's  properties  as  of  Fig. 5* 
D,,0  =  iOl  . 


D 

Fig.  10.  Drift  (horizontal)  velocity  of  a  cylinder  falling  towards 
an  horizontal  plane  for  various  initial  orientations  Q0 
(in  rods).  Same  cylinder's  and  liquid's  properties  as  of  Fig. 5 


10  1. 


11.  Angular  velocity  of  a  cylinder  falling  towards  an  horizontal 
plane  for  various  initial  orientations  (in  rads)  at  :D,C 

*  10  1,  5  1.  Same  cylinders'  and  liquids  properties  as  of  Fig. 5 


17-  Drift  (horizontal)  velocity  of  a  cylinder  falling  towards  an 
inclined  plane.  Same  parameters  and  conditions  as  of  Fig. 13. 


I  (»«<$ 

.  18.  Angular  velocity  of  a  cylinder  falling  towards  an  inclined 
plane.  Same  parameters  and  conditions  as  of  Fig.  12.  (D4 
reaches  a  minimum  at  the  point  P) . 


Fig.  19.  Angular  velocity  of  a  cylinder  falling  towards  an  inclined 
plane.  Same  parameters  and  conditions  as  of  Fig.  13-  ( 
reaches  a  minimum  at  the  point  P.). 

Motion  within  a  Flow 


kirc) 


Fig.  20.  Change  with  time  of  the  vertical  velocity,  A  (in  1/sec) .angular 
velocity  B  (in  rads/sec) ,  distance  from  the  wall  C  (in  1  )  ,  and 
orientation  anglelKin  rads)  of  a  cylinder  within  a  Couette  flow. 

X>,,*  =20  1. 

Same  cylinders  and  liquids  properties  as  of  Fig.  5. 


21.  The  vertical  velocity  A  (in  1/sec),  angular  velocity  B 

(in  rods/sec) ,  distance  from  the  wall  C  (ini),  and  orientation 

angle  D  (in  rods)  as  a  function  of  time.  Medium-air  of 
-1* 

«  1.718  X  10  poise;  particle  -  1  ■  20  j  b  *  0.5  p.->n  - 
$r  =2.23  g.  cm  channels  half  wiith  a  *  0.5  cm; 
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'*■  Discussion 

The  results  shown  in  the  above  figures  should  be  discussed  while  Keeping 
in  mind  the  behaviour  of  a  cylinder  under  creeping  motion  conditions  in  an 
unbounded  (far  from  walls)  fluid.  In  the  latter  case,  an  homogeneously  dense 
cylinder  which  settles  down  in  an  external  field  devoid  of  orienting  forces 
moves  with  a  neutral  stability,  viz.  a  fixed  orientation  angle  (12).  In  our 
case,  however,  the  cylinder  is  definitely  subjected  to  a  torque. 

Vertical  Wall;  Still  Medium 

Here  one  can  observe: 

i.  A  continuous  change  of  the  orientation  angle  6  caused  by  the  wall- 
created  torque, 

ii.  A  continuous  development  of  a  drift  and  an  increase  of  the  displa¬ 
cement 

It  turns  out  that  the  sign  of  the  torque  depends  on  the  angle  ©  ;for  every 
initial  distance  of  the  particle  from  the  wall  there  exists  a  critical  angle  q* 
such  that 


L  F  <  0  when  1 8  |  >  q  *  , 

L,  >  0  when  |  <  f 

and 

-  0  when  18  1=  © 

The  correspondence  between  theory  and  experiments  depicted  in  Figs.  5  through  7 
is  really  good,  which  produces  a  confidence  in  our  calculatlonal  method. 


y  v 
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Horizontal  Wall;  Still  Medium 

Here  particles  motion  is  always  two-dimensional.  Its  characteristics  are- 

i.  A  continuous  decrease  of  the  sedimentation  velocity,  ,  as  the 
plane  is  approached, 

ii.  An  enhanced  wall  influence  for  small  6.-  , 

i i i .  A  negl igigble  effect  for  large  c  >  0 .  ‘O  : 

iv.  An  increased  value  and  sensitivity  to  C 0  of  the  angular  velocity 
as  the  particle-plane  gap  is  reduced, 

v.  A  constant  drift  velocity,  ,  till  very  small  distances  from 

the  pi  >ne. 


It  is  interesting  to  compare  these  results  with  those  of  a  spherical  particle-  i 

horizontal  wall  situation(2)  where  a  similar  behaviour  is  observed. 


Inclined  Wall;  Still  Medium;  Two-dimensional  Situation 

From  the  results  for  the  inclined  wall  it  is  seen  that: 

\ 

i.  In  the  case  of  a  vertical  initial  orientation,  the  torque  is 
positive,  which  will  cause  the  cylinder  to  drift  from  right  to  the  left  (Fig. 12). 

ii.  In  the  case  of  an  horizontal  in;tial  orientation,  the  torque  is 
negative  but,  neverthel ess , the  cylinder  will  drift  in  the  former  direction  too 
(Fig.  13). 

iii.  In  the  case  of  a  wal 1  - incl i nat ion  of  a  little  less  than  75°  for  a 

vertical  initial  e*y) i nder-or i entat ion  and  of  a  little  less  than  60°  foi  an 

horizontal  wal 1  - i ncl i nat ion ,  the  angular  velocity  w  and  the  horizontal 

separation  02  respectively  reach  a  minimum  and  a  maximum  at  the  same  time. 

Thus,  the  particle  will  not  collide  with  the  surface,  which  is  of  a  great  signifi¬ 
cance. 
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Here  one  should  note  the  aerosoi  specificity  as  compared  to  the 

motion  within  our  liquid.  Since  both  the  mass  ar.d  the  moment  of  inertia  of 

2 

a  cylindrical  particle  are  proportional  to  b  ,  :t  comes  out  that  the  thinner 
is  the  particle  the  greater  are  its  translational  (in  l/sec)  and  rotational 
velocities  attained  after  a  certain  time.  However,  if  one  returns  to  dimen¬ 
sional  quantities,  one  observes  that  the  overal 1  particle  drifts  in  aerosol 
and  liquid  systems, as  caused  by  a  plan  ar  wall,  will  still  be  the  same.  For  an 
aerosol  system,  the  influence  of  the  wall  is  manifested  already  at  distances 
of  about  50  1  (  2)  . 

Motion  within  a  Flow;  Two-dimensional  Situation 

From  Fig.  20  it  can  be  seen  that  the  behaviour  of  our  cylinder  as  it  rotates 
in  the  gradient  flow  is  essentially  the  same  as  that  far  from  a  wall  ("Depo¬ 
sition  from  Flow1',  Final  Report  DAERO-76-G- 1 4 ,  1978,and  Gallfly  and  Eisner, 

J.  Colloid  Interface  Sci .  68^  320-37  (1979)).  The  wall  effect  is  shown  In  the 
decrease  of  the  second  maximum  of  the  vertical  velocity  V,  as  the  particle 
nears  the  horizontal  plane  (curve  A);  the  effect  is,  however,  over -shadowed  in 
the  angular  velocity  B  curve.  It  is  noteworthy  that  the  (rotating)  particle  spends 
more  than  85%  of  its  time  in  an  horizontal  orientation. 

?9l5euli le_Flow 

The  behaviour  in  this  case  is  essentially  the  same  as  that  within  a  Couette 

flow. 


i.  For  cylindrical  particles  there  is  a  pronounced  wall  effect  which 
may  be  of  a  great  importance  ir,  a  motion  within  small  cavities  like  those 
of  the  human  lung,  and  near  small  aerosol  collecting  obstacles. 

ii.  The  sedimentation  efficiency  of  the  particles  on  an  inclined  or 
horizontal  plane  is  decreased  due  to  the  effect. 

iii.  The  surface  concentration  of  the  particles’  deposit  on  an  inclined  plane 
has  a  greater  dispersion  as  compared  to  that  for  which  the  wall  does  not 

exert  any  influence. 

iv.  The  collision  of  the  particles  with  a  vertical  wall  is  very  small  due 

to  their  "repulsion"  from  it;  only  particles  with  O  <  6  will  have 

a  slight  chance  of  colliding. 

Finally,  extending  the  results  to  motion  near  curved  convex  surfaces,  one  can 
assume  that  the  smaller  Is  the  angle  between  the  tangent  plane  and  the 
surface,  the  greater  is  the  (still  medium)  collision  efficiency  which  is 
defined  as  the  ratio  between  the  number  of  settled  particles  to  the  number  that 
would  have  sedimented  had  all  trajectories  been  vertical. 

V.  Appl ications 

The  slender  body  method  of  calculating  the  fluiddynamic  force  and  torque 
on  the  particles  seems  to  be  more  useful  than  that  based  on  Oberbeck's  (4) 
and  Jeffery's  (5)  equations  since  it  does  not  depend  on  any  limiting 
assumptions  concerning  the  flow  in  which  the  particle  is  immersed.  For  example, 
any  numerical ly  (computer)  solved  field  of  flow  can  be  introduced  into  the 
slender  body  equations  where  It  is  necessary  to  know  the  undisturbed  velocity 


r 
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(Eq.  119]  )  at  the  middle  of  each  of  the  N  divided  cylinders  sections. 

This  has  been  performed  (with  no  wall  effect .however)  in  a  DFG  (Germany) 
supported  study  in  calculations  which  checked  quite  satisfactorily  with 
experiments  (13). 

B.  BROWNIAN  DIFFUSION 

I.  Theoretical 
1 .  Method 
a.  General 


The  diffusion  of  aerosol  particles  in  air  is  generally  controlled  by  a 
Brownian,  a  turbulent  or  a  mixed  wechan Ism;  near  surfaces,  where  turbulence 
subsides,  it  is  dictated  mainly  by  the  first. 

In  the  case  of  spherical  particles,  there  has  been  already  much  study  on 
the  numerous  aspects  of  Brownian  diffusion,  but,  very  little  has  been 
performed  with  regard  to  the  process  for  the  nonspherical  counterparts .  The 
difficulty  arises  from  the  general  non-isotroplc  nature  of  teh  fluiddynamic 
interaction  between  particles  and  medium  which  makes  it  necessary  to  consider 
the  spatial  orientation  and  the  rotation  of  the  former. 

A  useful  way  of  treating  the  problem  is  based  on  Brenner's  formal  structure  (14) 
in  which  one  assumes  a  six-dimensional  location-orientation  hyperspace  where 
the  general  number  concentration  Is  given  by  the  solution  of  the  conser¬ 
vation-equation 


w  *. 
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ensemble. 
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In  Brenners'  (formal)  treatment  (14)  and  for  an  axi -symmetric  particles'  case, 
one  has  instead  of  [27] 


— -  +  +  .  A 

9t  3r  w  ^  o 


O 


[30] 


where 


J-  ae 


C  u 


,and 


C(rte  ,*)  =  c  (r(i  )  f  Cir,e,-t) 


[31] 


In  that  case,  one  obtains  for  the  orientation  distribution  density  f  the 
governing  equation 


Trf/'di  +  ~Ti '  -  X  ( 


C32] 


which  takes,  for  example,  the  form 


~dC 

3+ 


f  r*. 


[33] 


where  fc», 0  are  the  polar  angles  and  the  rotational  Pecl&t  number  rPe 
is  defined  by 

H  =  4M-  •  0<a 


For  a  quiet  medium,  the  orientation  density  f  has  been  calculated  to  be 
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For  a  flowing  medium,  the  convective-rotation  equation  of  axi symmetric 
particles  for  example,  viz. 


Xf  -  V,  •  C  +  M ) 


[36: 


has  been  (firstly)  solved  by  Peterlin  (16)  when  the  Couette  conditions  and  the 
pertinent  Jeffery  equations  for  Qj  (5)  hold.  The  solution  shows  f  to  be  a  sole 
function  of  the  particle  Peclet  number  . 


b.  The  Aerosol  Specificity 

From  the  above  it  can  be  observed  that  the  (crucial)  density  function  f 

r 

for  axi-symmetric  particles  depends  on  Dx  and  to  .  The  latter  have  been 
calculated  in  the  case  of  a  liquid  suspension  by  classical  fluiddynamic 
methods  (17,5);  however,  in  the  case  of  aerosol  systems  this  Is  not  adequate. 

The  reason  for  the  inadequacy  stems  out  of  the  fact  that  in  the  very  instance 
of  the  only  (small)  aerosol  particles  which  show  a  pronounced  Brownian  motion, 

the  medium  cannot  be  considered  as  continuous  any  more.  Thus,  both  the  rotational 
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diffusion  coefficient  X>A  and  the  gradient-rotational  velocity  of  the 
particles  have  to  be  calculated  from  first  principles  by  statistical-mechanics 
means . 

c.  Immediate  Aim 

It  was  intended  to  study  (for  the  first  time)  the  rotational  diffusion 
coefficient  of  cylindrical  or  discoidal  aerosol  particles  which  form  a  good 
geometrical  model  for  many  systems.  The  task  has  been  accompl i shed ,as  reported 
now  and  also  published  under  the  title  “On  the  Stochastic  Nature  of  the  Notion 
of  Nonspherical  Aerosol  Particles  ill.  The  Rotational  Diffusion  Diadic  and 
Appl icat ions"by  Alfred  D.  Eisner  and  Isaiah  Gal  lily,  J.  Colloid  Interface 
Sci.  8^,  214-33  (1981). 

The  gradient  rotational  velocity  of  the  particles  in  a  simple  (though  not 
limiting)  shear  flow  has  been  calculated  too  and  presented  in  the  publication: 

“On  the  Stochastic  Notion  of  Nonspherical  Aerosol  Particles.  IV.  Gradient 
Convective  Rotational  Velocities  in  a  Simple  Shear  Flow"  by  Alfred  0.  Eisner 
and  Isaiah  Gal  lily,  J.  Colloid  Interface  Sci.  1982  (in  press),  which  is  not 
included  here. 

Finally,  a  third  (summarizing)  article  is  in  preparation. 

d.  Medium,  Particle,  and  General  Governing  Equations 

Out  of  the  generalized  Einstein  equation 

“m  *  V  i«T  (*=t,»*,C)  f37] 
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between  the  translational,  rotational  or  coupling  particle  mobility  and  the 
respective  diffusion  tensors,  and  from  the  relation  to  the  operating 
torque  exemplified  by 


tF  =  -  m  •  cj  ,  [38] 

one  notes  that  the  core  of  the  calculations  is  imbedded  in  the  evaluation  of 
L  F  . 

The  (considered)  particles,  assumed  to  be  much  smaller  than  the  mean  free  path 
l  of  the  gas  molecules  (  K-n  >>  i  )  .are  taken  to  translate 
with  a  linear  uniform  velocity  with  respect  to  an  Inertia)  system 

of  coordinates  CK,  (fig.  22),  and  rotate  with  a  uniform  velocity  with 

respect  to  the  center  of  mass  M  in  a  quiet  medium  devoid  of  external  torques. 
Thus,  the  source  of  Lr  is  supposed  to  be  only  the  stochastic  bombardment 
of  the  gas  molecules. 

For  our  medium  (and  structureless  molecules),  the  distribution  of  the  molecu¬ 
lar  velocities  is  Maxwellian,  i.e. 


N.t  dv?=  N  (h/r)*  C-Kv.”t  v®‘.  v“)J  dvS’rfW  f»J 

where  ( 1 )  ,  (2)  , (3)  respectively  stand  for  yx  f  2-j,  and  h *n/afcr'r. 

Concerning  the  inertial  system  04  itself,  there  Is  no  preferred  position  of 
its  placement  in  space  since  the  latter  is  assumed  (by  the  above)  to  be 
isotropic. 


-  40 


Fig.  22.  The  external  system  of 
coordinates  0_  . 


The  particles'  surface  condItion^  idopted  in  our  analysis  are  those 
of  Maxwell  (18)  according  to  whom  one  has  to  distinguish  between  a  diffusely 
reflected  fraction  <*  of  gas  molecules  with  a  Maxwellian  distribution 
conforming  to  the  surface  temperature  Ts  and  a  specularly  reflected 
fraction  (  f-a*  )  for  which 


-  2tl<  (H*  V,) 
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From  the  point  of  view  of  an  observer  at  Oj.  the  velocity  of  a  gas 
molecule  relative  to  the  body-point  0^(Ffg.  22)  is 


v,  =  •  coxf?  -  UM 


[41] 


Thus,  the  number  of  these  molecules  of  velocity  between  va  to  v.  *  iva 
hitting  an  elemental  surface  el  £  around  that  point  in  a  unit  of  time 
i  s 


•  exrC-Uvja3  ot\  ,  [42] 

the  momentum  transfer  at  that  time  is 

yr\  (  Vj  —  gj  y  R  —  uH  )  ^  d\  ol£  ,  C*3] 

and  the  torque  imparted  to  the  particles  in  the  It  direction 
i  s 

(  Lp-  i?)<  )  ~  -  Nm  (h/rr)3/2{(  [k.(Rx  (va  -  ox  R  -  u„))] 

S’  ** 

y  [CVj-wxw  -  uH)-n]  e*H|Vil  a’v*  o IS  [44] 


which  is  a  typical  governing  equation. 
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The  contribution  of  the  reflected  molecules  is  expressed  In  a  simitar 
fashion  but  with  the  exchange  of  h  and  N  of  the  velocity  distribution 
of  the  incoming  molecules  with  and  N  ^  of  the  diffusely  reflected 

ones  as  determined  by  the  requirement  that  (for  a  steady  situation): 


i.  The  net  flux  of  molecules  at  each  particles  surface  element  is  zero; 

ii.  The  net  total  energy  transport  at  that  element  is  zero,  too. 

The  total  torque  imparted  to  the  rotating-translating  particle  can  be  put  as 


)  =  0-«)  [(Lr-  le.  +  'J 


+ 


oc[(LP-k) 


McoO 


]  . 


[45] 


2.  Calculat ional  Procedure 
a.  A  Spherical  Particle 

This  has  been  chosen  as  a  test  case  for  the  proposed  method.  Taking  here 
the  "natural"  coordinates  system  and  (Pig-  23),  one  has 


R  =  r  (  si*>  e  cos  $  i  si-r,5  s.'-rtJ  +■  c Os6  k)  , 


Y!  =  Sin«  c-osf  i  •+  Sin  &  Si-rl  i  +  COS  6  h  > 


[47] 


and 


* 

— 

A 

V 

M 

v*/ 

[48] 


where  A,  the  transformation  matrix  is 


ft 


Fig.  23.  The  spherical  coordinate  system 


Due  to  the  space  isotropy,  it  is  possible  to  write 


W  =  W  k 


and 


V  i  +  to.  ^k) 
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Ci-»0 


A  sample  procedural  way  for  calculating  (  L  f  -  k  ) 
is  given  here  as  follows. 


i  .  The  integrand 


_  fc  .  [  q  x  (  Vj  -  OJ  x  R  -  V „)]  -  u  x  R  -  UM 


)*<n] 


=  r  sin  6  vi  f  U n  Sl 


oosfi  v/-  t»*«**A  v/ 


U„  cos  0  sin1# 

ii.  Integration  boundaries 


[52] 


cv,--«  6  o  and  so  , 


as  n  =  t 


( .  n  =  vly‘  =  v,  _  cx>s  e 


[53] 


which  gives 


Also  , 


4  Vj.*"  4  Uw  c oi  & 


4v/(  V,.*  &  «*> 


[5A] 

[55] 


Hi.  Final  equation 


(L  -k f  f  f  f 


«r  f  »  f  _Um  Col  6 
I. 


y  e 


-h(vxr  +  v/tv/’)  ( J  I  dLv^Av/din*  r*.-n«cl# 


[56] 


,  the  velocity  components  transformation  Jacobian  is  1 
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and 


where  B,  the  transformation  matrix,  is 


Placing  the  (^(external)  coordinate  system  (Fig.  24)  so  that  the  polar 
axis  of  the  particle  will  coincide  with  I?  ,  and  assuming  an  arbitrarily 
directed  angular  velocity  o>  =  w,  l  +  cOif  +  co*  fe 
one  can  carry  out  a  piecewise  calculation  of  the  total  torque  imparted  to 
the  envelope  and  each  base  of  the  particle.  By  order  of  magnitude  conside¬ 
rations  (APPENDIX  2),  it  can  be  shown  that  the  ordinary  translational 
velocity  of  an  aerosol  particle,  ,  Is  negligible  with  respect  to 

and  c»>  x  R  (Eq.  [44]  ).  Consequently,  It  becomes  possible  to 

express  the  integrand  of  [44]  for  the  molecules  which  hit  the 
particles  envelope,  for  example,  as 
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,  i 

I  t  =  l  -  [  R  X  (Vj-  CJX  r(v^  -  O  x  R)  .  nj  =  -  i  iil\4  v 

+  ("  COi  2*-  CO,  ya  sivt^P  -  co,  z^in*#  f  CJ2  f  *0 Otp  SrnP  +  COi^COtP  Sm  * 

4.  CJ,,  J?cos  *)  V,*  +  (  CJ,  5  H  Sin^  -  Wj  J*  sin#  cot*)  /*%  C-  W,2a5»ri*  CM* 
+  COi  i2COSa  ♦  )  Vj.^  +-  ?  S/-M  *  ***»-»«**  *V  -CO.^3*!^* 

+•  CO,1  2  cos1*  *m  p  -  cjj1?1?  C0ta*  $i"n  *  -  5  2  i,.„  *  +  2  UJ,  CO,  J3'* 

o»  *  s-tj1*  +  o,  co,  c«s  <*  +  w,  co3  5  Z2  oos  *  sr-n  4>  -  o>,  co3  j  I-2 

cos1*  .  [63] 


The  integration  limits  in  the  (governing  equation)  [44)  are  set  to  be 


y,  •*  =.  v,J  +  w.is'-r*  -  Z  cos  *  -  o  ,  Which  leads  to 


-  co  L  v/  6  CO,  2  cos  *  -  CO,  2  s  m  *  t 


and 


4  ^ 

-«®  V,  J  v,  c  » 


! 


F i na Illy,  in  this  case,  the  torque  can  be  written  in  the  form 

(lf  •  t  ^  =  -  Nvn  (  b/rr)^2  ff  f  f  X,  e'h  ^  +  v**+ 

S't  **  *♦  »* 

cfv/  d\J*  dv*d&  [64] 

where  the  velocity  limits  are  given  by  the  above,  and  &  &  =  f  ot  4 


(0  *  *  k  ITT  J  '  -J-”  *  4  ~T"  ) 


iSKMltt 
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For  the  case  of  which  almost  always  occurs  (APPENDIX  2),  and 

when  utilizing  certain  integral  values  (APPENDIX  2),  one  obtains  in  the 
typical  procedure  shown  here  that 


(  ■  /,T° 
(uF-  Wc  a 


N-m  '  h/Tr)'7*  '[-  - (-  CO, 


-  co 


,  *2cosa$  +  cot  coj  zm  *  -  coj*aCwa<>  !w*+  cjjSZCfl**,) 


+  (co, Vi  «****•*»♦  +  w.Vai***-"'*23*'17*  -w,Vi  «**«>*♦ 


+  00,  cDi^cos.*  +  co,coa  ?a2si-*»N>  +  cJjCOaS^Z  Oos<P  siVf  +  W.c^  J*a 


Cfistsfaf  ~cOa«s  Jj^COS*#) wiz  <*>*♦- cu ,Z-s*i"$)-t- W" •£)<(£>  [6?] 
where  0(...)  denotes  higher  powers  of  W;  which  can  be 

neglected  (APPENDIX  2). 

1 1 .  Results 

1.  The  Spherical  Particle 

The  total  torque  on  the  spherical  particle  for  u>  =  ^  k  and 
i s  found  to  be 


<«-r 


, (t«+) 

K)  =  -  ot  ( 


TT  \lfe- 

h  2 


CO 


N- 


(4 


8 

is 


K  ) 


[66] 


.a:  *v  w 
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which,  for  the  common  case  of  H  ^U  4X  1  (APPENDIX  2)  .turns 

to 


(L(r  .  k)^  =  <*  C4/»)(Wli)  ^  co  A/w  t-*  .  [67] 

Eg.  f67]  is  exact  ly  similar  to  that  obtained  by  Epstein  (l8)  in  a 
different  way  and  so  one  may  certainly  have  confidence  in  the  calculational 
method. 

It  should  be  remarked  that  in  the  above  case 


[68] 


as  can  be  easily  seen. 


For  the  arbitrary  UM  direction,  it  is  found  that 


(t 


F 


It  =  o<  (ir /  li  ]  W  N  yn  r1* 


xC-^-  +  _A_(i+  s.V.N')  K  UM4]  [69] 

which  becomes  identical  with  [66J  as  V  =  0  . 

2.  The  Cylindrical  or  Discoidal  Particle 
The  total  torque  on  these  particles  Is 
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(tr-  o£w,=-(t)'"c o,n™sl3  rc^-i  f  ♦  <»*♦  p’) 


a  C^i-  .ii.i)1*  tt-4  .s)  _  (  e 

^*g  iP  +  r  ;  +  ^co,  '  l* 


*  h*  o; 


Cu  >  CJj  /  £ 

8  M,  V  3 


- ra  )  +  *»  («•>*•*■  w») 


l  ^  _ r_ 

80  5*  L 


+  -J^tt)]  +  06-0  170] 


which  can  be  written  to  a  very  good  approximation  (APPENDIX  2)  as 


(l r •  i)"s  -<■£)''*«, 


+  V  *<8 


a  _  _L«V  JE=A_  a*' 

i,  r  ir+  s  r  > 


where 


Similarly  to  Eq.  71  ,it  is  possible  to  get  an  equation  for  the  torques 
imparted  to  the  rotating  aerosol  particle  along  its  other  axes;  these  come 
out  to  be 


(Lf  •  i  )  rv  -  U>: 


.  ct»-0 

(lf  •  W  )  ^  -  «s 


Thus,  the  (vector)  torque  can  be  written  as 


~r  _  [ 

where  J1  =  M  is  a  symmetric  diadic  whose  eigenvectors  lie 

along  the  symmetry  axes  of  the  particle.  A  similar  deduction  can  be 
reached  from  Landaus  and  Lifshitz's  (19)  reasoning  based  on  irrever¬ 
sible  thermodynamics  principles. 


a.  The  rotational  diffusion  coefficient 
From  Eq.  [38]  one  has  that 


a)  i 


(Lf. 


hT 


which  gives,  if  one  uses  the  kinetic  theory  relationship 


(75] 

| x  (*/") 


D.  = 


kT  K-n 


-V-P3)!! 


[76] 


whe  re 


Kn  =  T/f  . 


Eq.  [76]  should  be  compared  to  that  of  Edwards  (20) 
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as  is  presented  in  Fig.  25. 


Fig.  25.  The  rotational  diffusion  coefficient  of  Cylindrical 
particles  in  air.  ,Bj  ,Cj  -  values  of  this  study; 
ArB,C  -  Edwards'  values. 


III.  Discussion 

It  must  be  re-noted  here  that  our  calculations  apply  to  aerosol 
particles  whose  dimension,  in  any  direction,  Is  much  smaller  than  the 
mean-free  path  1  of  the  gas  molecules  (-0.069  f*w»  at  normal  conditions) 
and  to  a  quiet  medium.  For  larger  particles  within  a  flow,  one  has  to 
recur  to  the  Boltzmann's  integro-dlfferentlal  equation  of  the  velocity 
distribution  whose  (closed  form)  solution  Is  not  known  as  yet. 


The  results  of  Dx  themselves  show  that: 

i.  The  values  of  the  coefficient  are  of  the  order  of  10^-10^  sec 
which  makes  the  typical  relaxation  times  for  Brownian  rotation  phenomenta, 

-  =  :  /  i_  negligibly  small. 

ii.  The  calculated  values  are  greater  by  a  factor  of  about  30 

than  those  obtained  according  to  Edwards  in  classical  fluiddynamic  means,  which 
is  reflected  by  the  rotational  Peclet  number  assigned  to  a  diffusional  con¬ 
dition  of  orientable  particles  and,  hence,  by  the  orientational  distribution 
densi ty. 

IV.  Applications 

r 

As  an  example  for  the  use  of  in  a  real  problem,  it  is  possible 

to  take  the  diffusional  deposition  of  cylindrical  aerosol  particles  on  an 
infinite  plannar  surface  and  in  a  quiet  medium. 

a  General 

Employing  the  cartesian  coordinate  system  shown  in  Fig.  26,  one  can 
write  down  Eq.  [27] in  the  form 


||  =  *  ( V,- k)  Sf^e  oo ■££  +  c 


.+  +_ 


Si7.le  stm3*:  114-  +  CV  +  -Sf* 

o  x.  x  °  s 


si-*  4>  cos  — 5 — £ -  +  strip  oos6  SiVi# 

a  xz 


x 

+  S'n0cot0  cos*  1  +  5d  f-l—  JL_ ^  [78] 


-  5*  - 


Fig.  26.  The  cartesian  system  of  coordinates  used  in  the 
diffusion  analysis. 

where  zz  „  }  x  are  the  translational  diffusion  coefficients  of  the 
particle  parallel  and  perpendicular  to  its  polar  axis,  respectively.  Now, 
if  it  is  assumed  that  particles'  orientation  density  is  Independent  of 
thei r  location  viz. 


<T  (r,4>,  t)r  C  O.  ,  X,  ,  x3)  ♦  ,  O  , 

then,  instead  of  C78]  ,  one  obtains  the  two  equations: 


[79] 


whose  solution  is  given  In  Eq.  r35],and 


.  n 

TiTs  T$r  J 


ceo: 
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|f-  =  C  +  tC  ^4  (  1),,  -  J)j 


3X,1 


S i<ri*0  *in4$]  +  C  X> .  +  (  J3„  -*13  A  )  tOSXe  ] 

a*a 


d2  C 

¥*? 


+  i(-E>„-PA)[  »l7,2e  sin  <p  0OSi<fc  ^ — — 


5x,  dXj. 


Sni6  oo sff  sfr>  ^ 


T>zc 
3X*  »X3 


+ 


Sir©  c-»s©  Cos 


B^C 

0x3  3x, 


3 


[81] 


At  this  stage,  the  orientation  averaging  being  defined  by 


t r 


<3> 


^  ^A)a-ne  da  d  <P 
r  j  |  (o,  4> }  t )  S'lfl  o/«  <=/♦ 


one  reaches  the  final  equation 

jtt  ir 

c)C 


,  sc  v  S  S  “IW 

'  3t  '  - 


a[  -f  fa,!**  ,*)  seeded* 


■3c 


=  Cx>A  +  ('xi„  -Vl  )<*;«».*©  cov2*)]  ~fr  +  C  *x>tL+(3>u-  33  J 

o  x( 


<s,Va*.V4»>]|^-  +  [  j^  +  cVl+Cx*.  -^O<cosa«>]  [62] 

^  ‘'i  tf  X» 


since  the  mixed  terms  3AC  /  dx;  a  Xj  drop  out. 


b.  Specific 


In  the  specific,  above  diffusion  case  where  It  is  taken  that  concen¬ 
tration  gradients  exist  only  along  Xa 


3  c 
dt 


[83J 
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where  the  orientation-averaged  translational  diffusion  coefficient  ab  (t) 


i  s 


Jb  (t)  =  +  C  B.,  -  I.V4>  .  [84] 

Reforming  the  orientation-averaging  of  the  function  si™*©  ti-n4  ♦  one 
gets  that 


5b(t)=X  +  rfl3 

Concerning  the  solution  of  [83]  ,  it  is  advantageous  to  define  a  (new) 
vriable  P  by  (21) 

t 

P  =  /  «Z>  CfJcit’  [86J 

whereby  one  obtains  that  (for  c  =  c<>  at  t  3  0  and  c  *  o  at 

xa  =  o  and  t  >  o  ) 

c  =  c,  $  C  i  C  TrfJ^rr^1  * 


$  being  the  Error  Function. 

The  total  number  of  deposited  particles  M^after  time  t  will  be  then  (22) 
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M*  =  1  Vo  ( 


J  £  (V 


rr 


[88] 


c.  Results 


The  number  M  t  vs.  time  is  presented  in  Fig.  27  for  the  two 
cases : 

A.  When  and  +  DX  are  calculated  according  to  classical  continuum 

theory,  viz. 


t 

15 


H 


_ kr _ 

i  rr  ft  L  [  i/<fn  ( 1/  p;  - 


[83] 


and 


t 


_ k t _ 

h  ’iT  p  L  ['/<  C'/P)  +  C.S->] 


and  Edwards'  value  of 


r  E 
^  ) 


is  assumed.. 


[50] 


B.  When  these  coefficients  are  given  by  Dahnekes  molecular  regime  equations, viz. 


X>„ 


_ kT  K-y, 

Tr  f*  +  (i/p.)-  &-rr)/4>] 


[91] 


and 


•r 


_ kT  k'-v, _ _ 

*C(l  /0  +  «<l*n/0-  (i-TOA'/J 


[92] 


and  our  value  of  JD  ^ 


Is  taken. 
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Also,  M_  is  calcuaited  for  case  C  of  a  volumetrical  ly  equivalent 
spherical  particle. 

From  Fig.  27  it  is  seen  that  the  values  of  case  B  are  significantly  higher 
than  those  of  case  A  but  that  those  of  case  C  are  higher  still,  which  is 
resonable  if  one  remembers  that  diffusional  particle  mobility  is  a  property  of  the 
surface  of  the  latter  and  that  the  surface  of  the  volumetrical ly  equivalent 
sphere  is  the  smallest  of  all  forms. 

It  should  be  noted  that  in  the  st i 1 1 -medi um  diffusion  discussed  here 
is  essentially  independent  of  the  (very  large)  value  of  «£>  (*)  for  all  experimen¬ 
tal  times  of  interest.  The  value  of  B  j.  will  matter  in  difusion  in  flow 
where  the  rotational  Peclet  number  (  _  \c&\  /  *x)  )  does 

count . 


I 

Fig.  27 •  The  total  number  Mt  vs.  t. 

A,  case  A*,B ,  case  £  ;  C  ,  case  C  cited  above  . 

?  =  O.  oo'fS  f**r,  ,  L  =0.1  t  T  =  118.  (6  eK  ,  fi-  1f<*.  22  X  I0'i  poise 

I  =  6.  69  x  to"*  On  . 


C.  APPEND  1 X  1  ,2 


1 .  Append i x  1 

a .  Green  functions 

Generally,  expressing  the  disturbance  due  to  the  slender  body 
by  (Fig.  1) 


where 
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Cr')  -  <$  f‘4j  O',  f')  ol*  , 

f ^  (r1)  -  $  ?'  ?*•  fj'  (r1,  O  ol  *  tiz. 

s  i  nee 

*• "  =  r 1  •+  f '  (  r '  i  -p  )  . 

b.  Definitions 

C ,  -  C ,  <^<>s  fl  Cj  sen  0  ^  C2  =  -  C,  sin  6  •*•  C2  C«4  £>  f  C3  =  C3  ^ 

H„=  2  COS  fl  +  s.'oiafl  Mia  =  M2t  =  -  si-re  Cos  6  f  M,3=  =  0; 

M2i  =  2  s.n*0  -*•  cor^  ,  Mi3  *  Mw  -  o  }  M„  =  ±  . 

TT,i  —  K  co%l0  +  K22  sir-*©  j  =  “^^.za-tOsiVig  c/>»  6  } 

TT>i  ~  K23  s/™  ^  j  n~j.3  =  K  s.nJ0  •*•  *  a*  Oqs36  j  TT25  =  K2Ic<»<9, 

TTb|  =  Kw  s^e  ;  TTj2  =  K" OO-s  e  }  TTjj  as  I <gS  # 

^ii  ~  /  f-2)  Sin  6  ,  I  u  =  It(  =  C ?i  / f  *)  Sing  oos  6  ; 

1*»  =  X3.  =  (^?3/S3)  Sm  e  ,  x„  =  t*C0s*6  X33=  ($//  Ji)  . 

1,1  =  sinod  ,  a  SIH0  cosd  ,  r,,*  s  0  ,  *1,  *  001*6, 

=  Is'a  =  X,',  =  o. 

x„"  s*  oo*-*e  ,  X)a"  =  ia("  =  -sine  cose  ,  Xaa  2;»  =  x».-  ~ 


where  S  is  the  angle  between  x  and  x 


'<  s  = 


r  4-.  £«;  _ 


— •  dr' 


J  IS  -  r 


.v_;r,  -  J  i-  (r,  ■"  W- 


c.  Elaborate  procedure 


This  procedure  is  based  on  Eq.  [20].  Employing  the  transformat  ions 


x  - - >  X  y,  =  >7  A  ,  ,X  1 - »  y  r ,  =  •n’  A  whe.ee  fi  =  2  1/  b1  > 


one  defines 

F  1 )  =  A  (» -'/)')  ,  'ti')  =  ^  6  '-r  -  n’)  ,  m'!  =  0  , 

'■'*  Cti  ,  v»0  =  Si-nf  fn-n'/i  ,  V,  (  n,  77')  =  CXss  ©  ( -n -r1']  A  t 
r,( n,  y\')  =  0  ,  h,  (m,  m’)  _  sen  ©  C -n  -  77V  a  }  K  (tj,  -*■)  =  tj'A  *  ; 

<3  J  (  -n’  -  -n  -  0.S-)  A 

'i  °  (  n'  -  r  -  r.5)  A  T 
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and  writes  [17]  and  [)8]  in  terms  of  these  transformations. 

Thus,  u..  w  (x  x‘)  - >  ii  (>',71')  and  Eq.  [20]  is  expressed  as 

j  *  4  1 

u’ V r  |r 

r  6 

where  Jf'=  M  (  j- ■)  *•  *n‘  ,  f‘) 

=  iM(^r')^^.r)  ,  ^(r,r)= 


2.  Appendix  2 

a.  Order  of  magnitude  of 


Taking  a  typical  case  of  =  0(i)  and  calculating  h  to  be 

5".  ^7  x  io*  10  sec2  cm~z  for  T  =  2  98  .  1 6  *  K  and  ai  r  molecules, 

one  obtains  that 

h^L/M  =  0  ( *0  )  f  which  is  much  smaller 

than  unity. 

b.  Order  of  magnitude  of  U  v,  cot 

Taking  L  =  0.  I  x  io"*  c-m  ,  $  =  Q.OOlS  V  to'^crr}.^  lj. 

and  a  i  r  molecules  at  T  *  2*18  .  1 6  *  1C  ,^  =  18  4.22*  10“^  f>o  ise,  l  -6.6  » 


and  estimating  to,  by  the  law  of  equipartition  of  energy,  viz. 


-r  Ii  <«**>  =  t  *T  > 


rc.3] 


one  gets  that 


i>k 


D.  NOMENCLATURE 

a  -  constant  of  systems  geometry;  half-width  of  flow  channel 
A,B  -  transformation  matrices  defined  in  text 

b-  character i s t i c  equivalent  cross-section  radius  of  (slender)  particle 

b  -  local  equivalent  cross-section  radius 

r 

8.  ,  B."  ,  B''  ,  B.  -  constants 

c  -  particles'  total  number  concentration,  o  -  initial  value 

0  -  vector,  defined  in  text 
D  -  distance,  define.d  in  Fig.  2 

■f  <  «* 

D  i  0  .  D  -  translational,  rotational  and  coupling  diffusion  diadics, 
respect i ve ) v 

Oj.  -  rotational  diffusion  coefficient  around  a  mid-diameter 

V  j  -  translational  diffusion  diadic  components  parallel  and  perpendicular  to 
polar  axis,  respectively 
D.  , D  1  -  distances,  defined  in  Fig. A 

-  *0  :  +  (^D  ,  sVe  &:?•*<*> 

e  -  unit  vector  along  svmmetry  axis  of  particle 

f  -  orientation  distribution  density 

-  5tokeslets'  distribution  density;  f ^ j,  C )  -doublets’  distribution 

dens i ty  etc . 

-j-  (r"  -  drag  force  per  unit  area  of  particle 

v 

Fe  -  external  force;  F,  -  fluiddynamic  force 

g  -  acceleration  of  gravity 

h  -  particle-wall  distance;  also,  m/2kT 

I  ,  I;  -  moment  of  inertia  of  particle 
I  ,  I',  I"  -  diadics,  defined  in  text 
I,  -  interqrand,  defined  in  text 

^  -  rotational  flux  of  axi -symmetric  particles 

J  -  particles'  flux 

IJ|  -  transformation  Jacobian 

k  -  Boltzmanns  constant 

k,  ,  k^-  constants 

Kvt  -  Knudsen's  number 

1  -  semi-length  of  slender  particle 

1  -  mean  free-  path  of  gas  molecules 

L  -  length  of  cylindrical  particle 

I  -  avtornal  tnrnilA'  I  _  _  _  jf  1  IjJ  HHvnami  C.  tOrOUP. 


6i> 


m 

Mt 

M 

T1 

n ,  n 


y  * 


P,f 

r< 

r.r',' 

r 

9. 


aM  Cr* 
v  <r) 


vi;  v, 


,-tt) 

V  (r) 


,  »  x5 


-  mass  of  particle 

-  total  number  of  deposited  particles  in  time  t 

-  diadic,  defined  in  text 

-  mob i 1 i ty  diadic 

-  unit  vector  normal  to  dS 

-  running  indices 

-  number  of  molecules  bitting  ur»i  t  surface  per  unit  time 

-  density  of  Maxwell's  velocity  distribution 

-  average  number  of  gas  molecules 

-  defined  in  text 

-  rotational  Peclef  number 

-  particle  rotational  Peclet  number 

-  associated  Legendre  function 

-  typical  length  of  particle 

r*  -  distances  defined  in  Fig.  1  and  2 

-  location  radius  vector  ,  v-  _  rKd>ui  o+  sflitre 

-  radius  vector  from  geometric  center  of  particle  to  a  point  on  its 

surface. 

-  particles  surface;  els  -  surface  element 

-  time 

-  Temperature 

r'/  -  etc.-  Green  functions  associated  with 

-  translational  velocity  of  particle 

-  disturbance  velocity  of  slender  body 

-  velocities  of  gas  molecules  in  the  0X  system  and  relative  to  dS 

respectively;v'  -  velocity  of  reflected  molecules. 

-  translational  velocity  of  slender  body 

-  double-layer  potential 

-  position  distances  on  body's  axis  (Fig.l) 

,  X3  -  coordinates 


% 


Greek  letters 


-  traction  of  gas  molecules  reflected  diffusely 
i.  -  ?  w 

-  running  ino’ces 
i  -  Kronickers'  delta 

A  -  l/N 

<i-  Levi-Civita  permu  1  at  i  on  tensor 
^  -  fluid's  vort ic i ty 
f  -  anqle,  defined  in  text 
M -  polar  angles 
•A.  -  vector,  defined  in  text 
1  -  dynamic  viscosity  of  gas 
TT  -  diadic,  defined  in  text 
f  -  position  vector,  defined  in  Fig.  1 

-  gas  density,  cylinders  raaius 

-  particles'  density 

0  -  number  concentration  in  location-orientation  hyper-space; 

a  iso,  [  f%,  ,'2  1/  s.'J"' 

$  -  orientation  (pseudo) vector 

$  -  Error  Function 

N*  -  angle  between  UM  and  to 

cu  -  rotational  velocity  of  particles;  to,  ,  co2  ,  co3  -components 
SI  =  ”m*' 

Subsc  r i pts 
M  -  center  of  mass 
o  -  initial  value 
c  -  cylinder’s  envelope 
Superscr ipts 

(in)  -  hitting  molecules 

c*  -  t,r,c 

(r)  -  reflected  molecules 

(tot)  -  total 
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